Study of Correlated Motions to Detect the Conformational Transitions of the Intrinsically Disordered Sheep Prion Peptide.

Intrinsically disordered proteins (IDPs) are known for their random structural changes throughout their sequence based on the environment. The mechanism underlying these structural changes is difficult to explain. All biological processes are known to follow the direction through which they act. A study of the correlated motion can help to understand the direction of the change. Herein, we introduced the multivariate statistical analysis (MSA) technique to study the correlated motion of the peptide. The correlated motion of the sheep prion peptide was studied with the change in the temperature and solvent. These techniques helped to identify the contributing residual motions that helped to form the different secondary structures of the protein and also the triggering factors that drive these sorts of residual motions. The structural details match the experimentally reported data. It was found that the direction of the change of the secondary structure for this peptide shifted from the C-terminal to the N-terminal with an increase in the temperature. It was found that the involvement of the hydrophobic residues present at the C-terminal and the middle residues (residues 12-17) is responsible for forming a ß-sheet at the normal temperature. Hydration water was found to play an important role in this change. Insights gained from this study can be used to design strategies for desirable structural changes in the IDPs.

Designing Reaction Coordinate for Ion-Induced Pore-Assisted Mechanism of Halide Ions Permeation through Lipid Bilayer by Umbrella Sampling.

Ion permeation mechanism through lipid membranes helps to understand cellular processes. We propose new reaction coordinates that allow ions to permeate according to their water affinity and interaction with the hydrophilic layer. Simulations were done for three different halides (F-, Cl-, and I-) in two different lipid bilayers, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dinervonoyl-sn-glycero-3-phosphocholine (DNPC). It is found that the involvement of the water molecules decreases the free energy barrier. The ions were found to follow different pathways for permeation. Formation of proper pores required a collaboration effort of the hydration shell water molecules and the hydrophilic lipid layer, which was favored in the case of Cl- ions. The optimum charge density and good water affinity of Cl- with respect to F- and I- ions helped to form the pore. The effect was prominently seen in the case of DNPC membrane because of its higher hydrophobic thickness. The umbrella sampling results were compared with other methods such as the Markov state model (MSM) and well-tempered metadynamics (WT-metaD).

Vetting of New 2,5-Bis (2,2,2-trifluoroethoxy) Phenyl-Linked 1,3-Thiazolidine-4-one Derivatives as AURKA and VEGFR-2 Inhibitor Antiglioma Agents Assisted with In Vitro and In Silico Studies.

The bioactivity of 1,3-thiazolidin-4-one derivatives with a 2,5-bis (2,2,2-trifluoroethoxy) phenyl moiety was computationally developed and evaluated. All of the synthesized thiazolidin-4-one derivatives have their chemical structures characterized using a variety of methods, including nuclear magnetic resonance (NMR) (1H and 13C), high-resolution mass spectrometry (HRMS), and Fourier transform infrared (FTIR) radiation. A human glioblastoma cancer cell line (LN229) was used to investigate the purified derivatives' antiglioma cancer efficacy. By using the MTT, colony formation, and tunnel tests, respectively, the in vitro cytotoxic and apoptotic effects of these compounds were assessed. Thiazolidin-4-one derivatives 5b, 5c, and 5e were discovered to have the best efficacy against glioblastoma cells out of all of these compounds. The derivatives 5b, 5c, and 5e were determined to have respective IC50 values of 9.48, 12.16, and 6.43 g/mL. Computation results showed that the bioactivity evaluations of the compounds were quite significant. The bridging -NH group forms a hydrogen bond with Glu 260 of synthesized derivatives 5b, 5c, 5d, 5e, and 5h. The vast majority of freshly developed compounds obeyed Lipinski's rule of five, which is in line with the results that the ADMET model predicted. Additionally, molecular docking evaluation and molecular dynamics simulation investigations against the proteins AURKA and VEGFR-2 were conducted for the synthesized compounds to incorporate both in silico and in vitro data. The findings revealed that almost all of the compounds had considerable binding to AURKA and VEGFR-2 residues, with binding affinities ranging from -9.8 to -7.9 kcal/mol. Consequently, the results of the biological investigations and the docking scores demonstrated that thiazolidinone molecule 5e containing 4-chlorophenyl substituent may be considered as a potential moiety for glioblastoma cancer treatments.

Exploring the Barriers in the Aggregation of a Hexadecameric Human Prion Peptide through the Markov State Model.

The prefibrillar aggregation kinetics of prion peptides are still an enigma. In this perspective, we employ atomistic molecular dynamics (MD) simulations of the shortest human prion peptide (HPP) (127GYMLGS132) at various temperatures and peptide concentrations and apply the Markov state model to determine the various intermediates and lag phases. Our results reveal that the natural mechanism of prion peptide self-assembly in the aqueous phase is impeded by two significant kinetic barriers with oligomer sizes of 6-9 and 12-13 peptides, respectively. The first one is the aggregation of unstructured lower-order oligomers, and the second is fibril nucleation, which impedes the further growth of prion aggregates. Among these two activation barriers, the second one is found to be dominant irrespective of the increase in temperature and peptide concentration. These lag phases are captured in all three different force-field parameters, namely, GROMOS-54a7, AMBER-99SB-ILDN, and CHARMMS 36m, at different concentrations. The GROMOS-54a7 and AMBER-99SB-ILDN force fields showed a comparatively higher percentage of ß-sheet formation in the metastable aggregate that evolved during the aggregation process. In contrast, the CHARMM-36m force field showed mostly coil or turn conformations. The addition of a novel catecholamine derivative (naphthoquinone dopamine (NQDA)) arrests the aggregation process between the lag phases by increasing the activation barrier for the Lag1 and Lag2 phases in all of the force fields, which further validates the existence of these lag phases. The preferential binding of NQDA with the peptides increases the hydration of peptides and eventually disrupts the organized morphology of prefibrillar aggregates. It reduces the dimer dissociation energy by -24.34 kJ/mol.

Growth Reaction of Gold Nanorods in the Presence of Mutated Peptides and Amine-Modified Single-Stranded Nucleic Acids.

Conformation of biomolecules like DNA, peptides and amino acids play vital role during nanoparticle growth. Herein, we have experimentally explored the effect of different noncovalent interaction between a 5'-amine modified DNA sequence (NH2 -C6 H12 -5'-ACATCAGT-3', PMR) and arginine during the seed-mediated growth reaction of gold nanorods (GNRs). Amino acid-mediated growth reaction of GNRs results in a snowflake-like gold nanoarchitecture. However, in case of Arg, prior incubation of GNRs with PMR selectively produces sea urchin-like gold suprastructures, via strong hydrogen bonding and cation-π interaction between PMR and Arg. This distinctive structure formation strategy has been extended to study the structural modulation caused by two structurally close α-helical RRR (Ac-(AAAAR)3 A-NH2 ) peptide and the lysine mutated KKR (Ac-AAAAKAAAAKAAAARA-NH2 ) peptide with partial helix at the amino terminus. Simulation studies confirm that a greater number of hydrogen bonding and cation-π interaction between the Arg residues and PMR resulted in the gold sea urchin structure for RRR peptide against KKR peptide.

2,5-Bis(2,2,2-trifluoroethoxy)phenyl-tethered 1,3,4-Oxadiazoles Derivatives: Synthesis, In Silico Studies, and Biological Assessment as Potential Candidates for Anti-Cancer and Anti-Diabetic Agent.

In the present work, a series of new 1-{5-[2,5-bis(2,2,2-trifluoroethoxy)phenyl]-1,3,4-oxadiazol-3-acetyl-2-aryl-2H/methyl derivatives were synthesized through a multistep reaction sequence. The compounds were synthesized by the condensation of various aldehydes and acetophenones with the laboratory-synthesized acid hydrazide, which afforded the Schiff's bases. Cyclization of the Schiff bases yielded 1,3,4-oxadiazole derivatives. By spectral analysis, the structures of the newly synthesized compounds were elucidated, and further, their anti-cancer and anti-diabetic properties were investigated. To examine the dynamic behavior of the candidates at the binding site of the protein, molecular docking experiments on the synthesized compounds were performed, followed by a molecular dynamic simulation. ADMET (chemical absorption, distribution, metabolism, excretion, and toxicity) prediction revealed that most of the synthesized compounds follow Lipinski's rule of 5. The results were further correlated with biological studies. Using a cytotoxic assay, the newly synthesized 1,3,4-Oxadiazoles were screened for their in vitro cytotoxic efficacy against the LN229 Glioblastoma cell line. From the cytotoxic assay, the compounds 5b, 5d, and 5m were taken for colony formation assay and tunnel assay have shown significant cell apoptosis by damaging the DNA of cancer cells. The in vivo studies using a genetically modified diabetic model, Drosophila melanogaster, indicated that compounds 5d and 5f have better anti-diabetic activity among the different synthesized compounds. These compounds lowered the glucose levels significantly in the tested model.

Understanding the role of water on temperature-dependent structural modifications of SARS CoV-2 main protease binding sites.

Thermally stable and labile proteases are found in microorganisms. Protease mediates the cleavage of polyproteins in the virus replication and transcription process. 6 µs MD simulations were performed for monomer/dimer SARS CoV-2 main protease system in both SPC/E and mTIP3P water model to analyse the temperature-dependent behaviour of the protein. It is found that maximum conformational changes are observed at 348 K which is near the melting temperature. Network distribution of evolved conformations shows an increase in the number of communities with the rise in the temperature. The global conformation of the protein was found to be intact whereas a local conformational space evolved due to thermal fluctuations. The global conformational change in the free energy ΔΔG value for the monomer and the dimer between 278 K and 383 K is found to be 2.51 and 2.10 kJ/mol respectively. A detailed analysis was carried out on the effect of water on the temperature-dependent structural modifications of four binding pockets of SARS CoV-2 main protease namely, catalytic dyad, substrate-binding site, dimerization site and allosteric site. It is found that the water structure around the binding sites is altered with temperature. The water around the dimer sites is more ordered than the monomer sites regardless of the rise in temperature due to structural rigidity. The energy expense of binding the small molecules at substrate binding is less compared to the allosteric site. The water-water hydrogen bond lifetime is found to be more near the cavity of His41. Also, it is observed that mTIP3P water molecules have a similar effect to that of SPC/E water molecules on the main protease.

Exploring the multiple conformational states of RNA genome through interhelical dynamics and network analysis.

The structural variation of RNA is often very transient and can be easily missed in experiments. Molecular dynamics simulation studies along with network analysis can be an effective tool to identify prominent conformations of such dynamic biomolecular systems. Here we describe a method to effectively sample different RNA conformations at six different temperatures based on the changes in the interhelical orientations. This method gives the information about prominent states of the RNA as well as the probability of the existence of different conformations and their interconnections during the process of evolution. In the case of the SARS-CoV-2 genome, the change of prominent structures was found to be faster at 333 K as compared to higher temperatures due to the formation of the non-native base pairs. ΔΔG calculated between 288 K and 363 K are found to be 10.31 kcal/mol (88 nt) considering the contribution from the multiple states of the RNA which agrees well with the experimentally reported denaturation energy for E. coli α mRNA pseudoknot (∼16 kcal/mol, 112 nt) determined by calorimetry/UV hyperchromicity and human telomerase RNA telomerase (4.5-6.6 kcal/mol, 54 nt) determined by FRET analysis.

Influence of Ion Specificity and Concentration on the Conformational Transition of Intrinsically Disordered Sheep Prion Peptide.

The structural sensitivity of the intrinsically disordered proteins with the ions has been observed experimentally; however, it is still unclear how the presence of different metal ions affects structural stability. We performed an atomistic molecular dynamics simulation of sheep prion peptide (142-167) in presence of different monovalent, divalent ions at various concentrations to find out the effect of the size, charge, and ionic concentration on the structure of the peptide. It is found that Li+ ions have a higher survival probability compared to Na+ , K+, and Mg2+ affecting the solvation structure of the protein leading to the alpha-helix structure. At high concentration, due to the increase in the ion-solvent and counter-ion interactions, the effect of the ions is screened on the surface of the protein and hence no ion specificity is observed. This study demonstrates how ions can be used to regulate the protein structure and function that can help in designing drugs.

Molecular mechanism of inhibition of COVID-19 main protease by ß-adrenoceptor agonists and adenosine deaminase inhibitors using in silico methods.

Novel coronavirus (COVID-19) responsible for viral pneumonia which emerged in late 2019 has badly affected the world. No clinically proven drugs are available yet as the targeted therapeutic agents for the treatment of this disease. The viral main protease which helps in replication and transcription inside the host can be an effective drug target. In the present study, we aimed to discover the potential of ß-adrenoceptor agonists and adenosine deaminase inhibitors which are used in asthma and cancer/inflammatory disorders, respectively, as repurposing drugs against protease inhibitor by ligand-based and structure-based virtual screening using COVID-19 protease-N3 complex. The AARRR pharmacophore model was used to screen a set of 22,621 molecules to obtain hits, which were subjected to high-throughput virtual screening. Extra precision docking identified four top-scored molecules such as +/--fenoterol, FR236913 and FR230513 with lower binding energy from both categories. Docking identified three major hydrogen bonds with Gly143, Glu166 and Gln189 residues. 100 ns MD simulation was performed for four top-scored molecules to analyze the stability, molecular mechanism and energy requirements. MM/PBSA energy calculation suggested that van der Waals and electrostatic energy components are the main reasons for the stability of complexes. Water-mediated hydrogen bonds between protein-ligand and flexibility of the ligand are found to be responsible for providing extra stability to the complexes. The insights gained from this combinatorial approach can be used to design more potent and bio-available protease inhibitors against novel coronavirus.Communicated by Ramaswamy H. Sarma.

Deciphering the competitive inhibition of dihydropteroate synthase by 8 marcaptoguanine analogs: enhanced potency in phenylsulfonyl fragments.

The emergence of sulfa-drug resistance and reduced efficacy of pterin-based analogs towards Dihydropteroate synthase (DHPS) inhibition dictate a pressing need of developing novel antimicrobial agents for immune-compromised patients. Recently, a series of 8-Marcaptoguanin (8-MG) derivatives synthesized for 6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (experimental KD ∼ 100-.0.36) showed remarkable homology with the pteroic-acid and serve as a template for product antagonism in DHPS. The present work integrates ligand-based drug discovery techniques with structure-based docking, enhanced MD simulation, and MM/PBSA techniques to demonstrate the essential features of 8-MG analogs which make it a potent inhibitor for DHPS. The delicate balance in hydrophilic, hydrophobic substitutions on the 8-MG core is the crucial signature for DHPS inhibition. It is found that the dynamic interactions of active compounds are mainly dominated by consistent hydrogen bonding network with Asp 96, Asn 115, Asp 185, Ser 222, Arg 255 and π-π stacking, π-cation interactions with Phe 190, Lys 221. Further, two new 8-MG compounds containing N-phenylacetamide (compound S1, ΔGbind-eff = -62.03 kJ/mol) and phenylsulfonyl (compound S3, ΔGbind-eff = -71.29 kJ/mol) fragments were found to be the most potent inhibitor of DHPS, which stabilize the flexible pABA binding loop, thereby increasing their binding affinity. MM/PBSA calculation shows electrostatic energy contribution to be the principal component in stabilizing the inhibitors in the binding pocket. This fact is further confirmed by the higher energy barrier obtained in umbrella sampling for this class of inhibitors.

Temperature-Dependent Conformational Evolution of SARS CoV-2 RNA Genome Using Network Analysis.

Understanding the dynamics of the SARS CoV-2 RNA genome and its dependence on temperature is necessary to fight the current COVID-19 crisis. Computationally, the handling of large data is a major challenge in the elucidation of the structures of RNA. This work presents network analysis as an important tool to see the conformational evolution and the most dominant structures of the RNA genome at six different temperatures. It effectively distinguished different communities of RNA having structural variation. It is found that at higher temperatures (348 K and above), 80% of the RNA structure is destroyed in both the SPC/E and mTIP3P water models. The thermal denaturation free energy change ΔΔG value calculated for the long-lived structure at higher temperatures of 348 and 363 K ranges from 2.58 to 2.78 kcal/mol for the SPC/E water model, which agrees well with the experimentally reported thermal denaturation free energy range of 2.874 kcal/mol of SARS CoV-NP at normal pH. At higher temperatures, the stability of RNA conformation is found to be due to the existence of non-native base pairs in the SPC/E water model.

Theoretical insights into molecular mechanism and energy criteria of PARP-2 enzyme inhibition by benzimidazole analogues.

The emergence of poly (ADP-ribose) polymerase (PARP) inhibitors targeting a class of PARP enzymes has gained a great interest in cancer therapy. Majority of the PARP inhibitors are not isoform-selective which may cause unwanted off-target effects. In the present study, we explore the molecular mechanism and energy requirements for PARP-2 inhibition. This involves docking studies, frontier molecular orbital analysis, 500 ns molecular dynamics simulation (MD), binding free energy analysis and principal component analysis. The results clearly suggest the importance of hydrogen bonding (Gly429, Gln332, Ser470, Tyr455) and π-π stacking interactions (His428, Tyr455, Tyr462, Phe463, Tyr473) between residues and the inhibitor. Presence of lowest unoccupied molecular orbitals favors π-π stacking interactions and highest occupied molecular orbital orbital favors hydrogen-bonding interactions in the ligands. The stability of most active/PARP-2 complex is confirmed by hydrogen bonding and π-π stacking interaction parameters. Molecular-mechanics Poisson-Boltzmann surface area energy calculations showed that van der Waals and nonpolar solvation energy terms are crucial components for the stable binding of the ligands. Per residue analysis showed that tyrosine, histidine, and phenyl alanine residues are responsible for hydrophobic interactions with the ligands. Four new inhibitors are designed based on this study and the stability of PARP-2/inhibitor complex is validated by MD, density functional theory studies, and ADME/toxicity properties. Information from the present study can serve as a basis for designing new isoform-selective PARP-2 inhibitors.

Epitope-Based Potential Vaccine Candidate for Humoral and Cell-Mediated Immunity to Combat Severe Acute Respiratory Syndrome Coronavirus 2 Pandemic.

The emergence of severe acute respiratory syndrome from novel Coronavirus (SARS-CoV-2) has put an immense pressure worldwide where vaccination is believed to be an efficient way for developing hard immunity. Herein, we employ immunoinformatic tools to identify B-cell, T-cell epitopes associated with the spike protein of SARS-CoV-2, which is important for genome release. The results showed that the highly immunogenic epitopes located at the stalk part are mostly conserved compared to the receptor binding domain (RDB). Further, two vaccine candidates were computationally modeled from the linear B-cell, T-cell epitopes. Molecular docking reveals the crucial interactions of the vaccines with immune-receptors, and their stability is assessed by MD simulation studies. The chimeric vaccines showed remarkable binding affinity toward the immune cell receptors computed by the MM/PBSA method. van der Waals and electrostatic interactions are found to be the dominant factors for the stability of the complexes. The molecular-level interaction obtained from this study may provide deeper insight into the process of vaccine development against the pandemic of COVID-19.

Structural and Thermophysical Anomalies of Liquid Water: A Tale of Molecules in the Instantaneous Low- and High-Density Regions.

Water is believed to be a heterogeneous liquid comprising multiple density regions that arise because of the presence of interstitial molecules and can be differentiated by their structure as well as the existence of hydrogen-bonded pairs with varying strengths. First-principles molecular dynamics studies were performed at six different temperatures to investigate the effect of temperature on the thermophysical, structure, dynamics, and vibrational spectral properties of the water molecules using dispersion-corrected density functional theory. The variation of properties like density, cohesive energy, and compressibility with a change in temperature produces a trend that matches with the experiments and resembles the experimentally observed anomalous behavior. We explore the possibility of explaining the trends in calculated properties by analyzing the structure and dynamics of the water molecules in terms of instantaneous low- and instantaneous high-density regions that are found during the simulation time. The dynamics of these two types of water molecules were studied by calculating the lifetime from the proposed autocorrelation functions. The lifetime of formation of instantaneous low-density water is found to decrease with an increase in temperature, whereas the lifetime of instantaneous high-density water is found to be maximum at 298 K among all the considered temperatures. The presence of more interstitial water molecules is observed at this temperature. The signature of these water molecules is found in the radial distribution function, spatial distribution function, void distribution, configurational space, orientational dynamics, and spectral diffusion calculations. It is also found that around 298 K, these water molecules are present distinctively that mix up with the first and second solvation shells with the rise of the temperature. The outlook of the reported results can be extended to other thermodynamic conditions to explain some of the anomalous properties, which can be related to the presence of the interstitial molecules in water.

Effect of hydrophobic and hydrogen bonding interactions on the potency of ß-alanine analogs of G-protein coupled glucagon receptor inhibitors.

G-protein coupled glucagon receptors (GCGRs) play an important role in glucose homeostasis and pathophysiology of Type-II Diabetes Mellitus (T2DM). The allosteric pocket located at the trans-membrane domain of GCGR consists of hydrophobic (TM5) and hydrophilic (TM7) units. Hydrophobic interactions with the amino acid residues present at TM5, found to facilitate the favorable orientation of antagonist at GCGR allosteric pocket. A statistically robust and highly predictive 3D-QSAR model was developed using 58 ß-alanine based GCGR antagonists with significant variation in structure and potency profile. The correlation coefficient (R2 ) and cross-validation coefficient (Q2 ) of the developed model were found to be 0.9981 and 0.8253, respectively at the PLS factor of 8. The analysis of the favorable and unfavorable contribution of different structural features on the glucagon receptor antagonists was done by 3D-QSAR contour plots. Hydrophobic and hydrogen bonding interactions are found to be main dominating non-bonding interactions in docking studies. Presence of highest occupied molecular orbital (HOMO) in the polar part and lowest unoccupied molecular orbital (LUMO) in the hydrophobic part of antagonists leads to favorable protein-ligand interactions. Molecular mechanics/generalized born surface area (MM/GBSA) calculations showed that van der Waals and nonpolar solvation energy terms are crucial components for thermodynamically stable binding of the inhibitors. The binding free energy of highly potent compound was found to be -63.475 kcal/mol; whereas the least active compound exhibited binding energy of -41.097 kcal/mol. Further, five 100 ns molecular dynamics simulation (MD) simulations were done to confirm the stability of the inhibitor-receptor complex. Outcomes of the present study can serve as the basis for designing improved GCGR antagonists.

Computational insights into factor affecting the potency of diaryl sulfone analogs as Escherichia coli dihydropteroate synthase inhibitors.

Dihydropteroate synthase (DHPS) is an alluring target for designing novel drug candidates to prevent infections caused by pathogenic Escherichia coli strains. Diaryl Sulfone (SO) compounds are found to inhibit DHPS competitively with respect to the substrate pABA (p-aminobenzoate). The extra aromatic ring of diaryl sulfone compounds found to stabilize them in highly flexible pABA binding loops. In this present study, a statistically significant 3D-QSAR model was developed using a data set of diaryl sulfone compounds. The favourable and unfavourable contributions of substitutions in sulfone compounds were illustrated by contour plot obtained from the developed 3D-QSAR model. Molecular docking calculations were performed to investigate the putative binding mode of diaryl sulfone compounds at the catalytic pocket. DFT calculations were carried out using SCF approach, B3LYP- 6-31 G (d) basis set to compute the HOMO, LUMO energies and their respective location at pABA binding pocket. Further, the developed model was validated by FEP (Free Energy Perturbation) calculations. The calculated relative free energy of binding between the highly potent and less potent sulfone compound was found to be -3.78 kcal/ mol which is comparable to the experimental value of -5.85 kcal/mol. A 10 ns molecular dynamics simulation of inhibitor and DHPS confirmed its stability at pABA catalytic site. Outcomes of the present work provide deeper insight in designing novel drug candidates for pathogenic Escherichia coli strains.

Interstitial Voids and Resultant Density of Liquid Water: A First-Principles Molecular Dynamics Study.

Many anomalous properties of water can be explained on the basis of the coexistence of more than one density states: high-density water (HDW) and low-density water (LDW). We investigated these two phases of water molecules through first-principles molecular dynamics simulations using density functional theory (DFT) in conjunction with various van der Waals-corrected exchange and correlation functionals. Different density regions were found to exist due to the difference in short-range and long-range forces present in DFT potentials. These density regions were identified and analyzed on the basis of the distribution of molecules and voids present. We defined a local structure index to distinguish and find the probability of occurrence of these different states. HDW and LDW arise due to the presence of "interstitial water" molecules in between the first and second coordination shells. The population of interstitial water molecules is found to affect the overall dynamics of the system as they change the hydrogen bond pattern.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins.

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

How Crystals Nucleate and Grow in Aqueous NaCl Solution.

Large-scale molecular dynamics simulations (64 000 particles) are used to examine the microscopic mechanism of crystal nucleation and growth in a slightly supersaturated solution of NaCl in water at 300 K and 1 atm. Early-stage nucleation is observed, and the growth of a single crystal is followed for ∼140 ns. It is shown that the nucleation and growth process is better described by Ostwald's rule of stages than by classical nucleation theory. Crystal nucleation originates in a region where the local salt concentration exceeds that of the bulk solution. The early-stage nucleus is a loosely ordered arrangement of ions that retains a significant amount of water. The residual water is slowly removed as the crystal grows and evolves toward its stable anhydrous state.